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Abstract 

The Melnikov method is apphed to periodically perturbed open systems modeled by an inverse- 
square-law attraction center plus a quadrupolelike term. A compactification approach that regular- 
izes periodic orbits at infinity is introduced. The (modified) Smale-Birkhoff homoclinic theorem is 
used to study transversal homoclinic intersections. A larger class of open systems with degenerated 
(nonhyperbolic) unstable periodic orbits after regularization is also briefiy considered. 
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I. INTRODUCTION 



Since the pioneering work of Poincare [0] in celestial mechanics in which the mathematical 
basis of deterministic chaos in compact phase space systems was laid down, the study of 
homoclinic phenomena in closed systems with hyperbolic unstable periodic orbits has allowed 
the understanding of a rich variety of nonlinear effects in physics, chemistry, and biology 
0. Due to its universality, models in which unstable periodic orbits are subjected to small 
periodic perturbations has become one of the main paradigms of deterministic chaos [^]. An 
analytical tool to study such models is the Melnikov method |^, ^ in connection with 
Smale-Birkhoff homoclinic theorem 0, ^ , and Kolmogorov-Arnold-Moser (KAM) theory in 
the Hamiltonian case |^. 

The Melnikov function describes the transversal distance between the stable and unstable 
manifolds associated to an unstable periodic orbit. Its isolated odd zeros indicate transver- 
sal intersections between these manifolds, and hence the onset of chaos [|10]- Examples of 
applications of the Melnikov method in gravitation are the motion of particles in perturbed 



two- and three-dimensional Stackel potentials [0, |T2[, the chaotic evolution of cosmologi- 



cal models |T^, the study of orbits around a black hole perturbed by either gravitational 
radiation |jT^ or an external quadrupolar shell |]15|, and the bounded motion of particles in 
a periodically perturbed attractive center described by a monopole plus a quadrupolelike 



potential were considered in Ref. |T^ . 

The Melnikov method has also been used in many other branches of physics. We find 
examples of applications of this method to the study of Josephson junctions |T^, |TB|, pla- 
nar periodic vortical flows [0, solitons [^, liquid crystals |2l[], and transfer dynamics of 
quasiparticles . 

Even for Hamiltonian systems fundamental questions about chaos in non-compact phase 
space systems remain to be answered. Among the more important unsolved questions are 
the notion of chaos itself and the lack of an adequate theory to deal with it. Partial results 



obtained in this area are the fractal techniques in scattering processes p3[, these are nu- 
merical techniques that present some difficulties due to the existence of different time scales 
for nearly bounded scattering. They are inadequate to the study of chaotic behavior arising 
from separatrices between bounded and unbounded orbits. Furthermore, they are unable 
to present a complete description of the chaotic motion as the one provided by analytical 
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methods in the closed system case. 

The aim of this paper is to study the homochnic phenomenon for a class of open systems 
that by a suitable change of coordinates can be approached in terms of an adequate for- 
mulation of the Melnikov method and Smale-Birkhoff homoclinic theorem. The change of 
coordinates regularizes the unstable periodic orbit at infinity and it compactifies the region 
of interest of the phase space; however, the phase space as a whole remains noncompact. 
Alas the resulting unstable periodic orbit is typically nonhyperbolic and the standard stable 



manifold theorem, needed to state the Melnikov method ||T0[, does not apply. McGehee ||24 
extended this theorem to degenerated cases in the context of the Newtonian three body 



problem. Xia and Dankowicz and Holmes [^, among others, used McGehee's result 
in connection with Melnikov method and Smale-Birkhoff homoclinic theorem to study the 
nonintegrability of the three body problem. 

Here we consider the equatorial motion of a particle moving in a potential described by a 
monopolar term plus a quadrupolelike contribution. This potential models the gravitational 
attraction of a galaxy bulge or any nonspherical celestial body; it also arises in general 
relativity in the study of the motion of a test particle around a Schwarzschild black hole, 
the quadrupole term being a general relativistic effect associated to the angular momentum 
of the particle in the reduced two-dimensional phase space, see for instance Ref. P7[] . 

In Sec. II the fixed saddle points associated to the monopole plus quadrupole system, as 
well as the coordinate transformation that regularizes these points at infinity, are studied. 
In the next section we present some mathematical preliminaries and the Melnikov method. 
The equations of motion are used to reduce this method to the analysis of simple graphics. 
We find that the perturbation induces transverse homoclinic orbits in some ranges of the 
parameters, and we apply the modified Smale-Birkhoff homoclinic theorem to verify the 
presence of a symbolic dynamics equivalent to a Smale horseshoe map; see Ref. In 
Sec. IV the study of the motion is completed with a presentation of Poincare sections that 
reveal different levels of chaotic behavior as a function of the parameters. Finally, in the 
last section, we make some remarks about the class of system in which the same kind of 
analysis can be performed. 
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II. THE HOMOCLINIC ORBIT AND PERTURBATIONS 



We shall consider the orbit of a particle in a plane under the influence of a force modeled 
by a potential with inverse square law plus a quadrupolelike term. It is convenient to work 
with dimensionless quantities. The motion of the particle is described by 

where r, p, t, Hq, and /3 are dimensionless quantities proportional to, respectively, the radius, 
the radial momentum, the time, the Hamiltonian function, and the quadrupole moment of 
the attraction center. The effective potential 

^e// = ^ ~ ~ ^ ^ ' (2) 

is presented in Fig. 1 for different values of the parameter (3 > 0. 

The natural space to study a periodically perturbed planar system is x S^, where 
the unstable periodic orbits have a proper meaning |jlO[. In the corresponding unperturbed 
autonomous case the phase space is in and the unstable periodic orbits reduce to fixed 
saddle points. The unstable periodic orbits are also reduced to fixed saddle points for the 
maps defined on Poincare sections M'^ x {to} C iR^ x S^. 

The above system presents a homoclinic loop associated to the hyperbolic fixed saddle 
point at (r,p) = (rA/,0), where 1/rM = 1/6/5 + ^(1/6/3)^ - 1/3/5, for (3 limited by 1/16 < 
P < 1/12. This case is important in the study of bounded orbits and was explored in Ref. 

In the present work we study the instabilities of unbounded orbits, the relevant values 
of the parameter are (3 = 1/16 and < f3 < 1/16. Let Hq = 0, if (3 = 1/16 the points 
r = tm = 1/4 and r = oo on the r axis represent fixed points such that the particle takes 
infinity time to reach to or to depart from each of these points. For zero energy still, if 
< /? < 1/16 the motion of the particle is restricted to the region between r = r_ and 
r = oo, where l/r_ = 1/4/3 — a/ (1/4/3)^ — 1/(3, and only the last point represents a fixed 
point in this range of values of (3. 

Since the orbits of interest are in a semi-infinity region bounded away from the origin, 
we can compactify this part of the phase space with a change of the position coordinate 
like r = with a < 0. We find that the transformation r = l/u"^ allows us to model 
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the problem in a way similar to Refs. p^, p5|, p6|. This new coordinate regularizes the 



fixed point at infinity that now is at the point {u,p) = (0,0). The zero energy orbits 
generate a heteroclinic loop for /3 = 1/16 associated to the hyperbohc fixed saddle point at 
{u,p) = (2,0) and the degenerated fixed saddle point at {u,p) = (0,0), and a homoclinic 
loop for < /5 < 1/16 associated to the degenerated fixed saddle point at {u,p) = (0,0). 
Degenerated in the sense that both eigenvalue of the linearized vector field are zero, as can 
be seen from the Hamiltonian equations 

du 1 

^ = -u* + u'-3fJu'. (4) 

The points {u,p) = (2, 0) and {u,p) = (0, 0) correspond to, respectively, the hyperbolic fixed 
saddle point at (r,p) = (1/4,0) for P = 1/16 and the degenerated fixed saddle point at 
(r, p) = (oo,0) for < j3 < 1/16. The homoclinic and heteroclinic loops are defined by the 
intersection between the stable and unstable manifolds on the {u,p) plane. The homoclinic 
loop for f3 = 1/18 is shown in Fig. 2 and the heteroclinic one is presented in Fig. 3. 

The explicit integration of the homoclinic and heteroclinic loops will be necessary to 
apply the Melnikov method and can be obtained from the first integral of motion {Hq = 0). 
For /3 = we find 

where v = 1/r, the time origin is take in the symmetry point of loop, and the sign refers to 
the upper (+) and lower (— ) parts of the loop. Analogously, for < /? < 1/16 it reads 

_ ^ /'2 / {2v+ + v.)Fi6, q) - 2{v+ + v.)E{5, q) 



3v^vT 



with 




arcsin^/^^^^)^^ ^ = A/ ~ ' (7) 

where F{S, q) and E{5, q) are elliptic integrals of first and second type in the Legendre 
normal form (Ref. [0, p. 224), and v± = 1/4/3 ± ^/{T/4py^^l/]3 are the roots of Veff = 0. 



Moreover, for (3 = 1/16 we get 

V2 f 1 1 \ 



t{v) = ± 



(t;3/2 (4/3)3/2) ^ 2v^ ^4/3 



8^2 



' m 



2 + v^y \2-^/I/3 J 



(8) 



with the choice t(4/3) = 0, where v = 4/3 is the local minimum of K//- 

Now, let us consider the Hamiltonian (|T|) perturbed by a periodic multipolar term of the 
form 

H = Ho + eHi, (9) 

Hi = r-"" cos{nt) (n>2), (10) 

where n = 2 is dipolar, n = 3 is quadrupolar, etc. These perturbations can model the 
attraction due to a distribution of masses with periodic motions that are placed inside the 
planet orbit. 

We shall consider our attraction center with a fix total mass. In other words, we excluded 
the monopolar case {n = 1) that represents a periodic variation of the mass. 

In the next sections we study how these perturbations can affect the dynamics of the 
system. 

III. MELNIKOV METHOD 

Powerful tools to study near integrable systems are Melnikov type of techniques that 
detect transversal intersections between the stable and unstable manifolds associated to a 
unstable periodic orbit. The presence of such transversal intersections is a guarantee of 
complicated dynamics and in some cases leads to a symbolic dynamics equivalent to the 
Smale horseshoe |]10| . 



In order to simplify the analysis it is convenient to abstract a little from the particular 
problem presented above. In what follows the loops are on the X plane and the manifolds 
are in the {X,6) space, where 6 = t mod 2tt/Q so that {X,6) G x S^, and denotes 
the section 6 = Oq. We consider a Hamiltonian of the form 

H{X,t) = HoiX)+eHi{X,t) (ifi 27r/f]- periodic in t), (11) 
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where Hq is integrable with homochnic (heterochnic) loop F associated to some hyperbolic 
fixed saddle point (s). Under hyperbolicity hypothesis it can be shown that for sufficiently 
small e the invariant manifolds are only deformed, and possibly their intersections become 
transversal, see for instance Ref. [|10| . 



Let Xq be a point on T and X*/"(6'o,£) be points on the stable/unstable manifolds such 
that they are on S^g, in the line perpendicular to F x {6q} at Xq and whose trajectories take 
the least amount of time to reach/depart any small neighborhood of the unstable periodic 
orbit. A computable measure of the transversal distance between the stable and unstable 
manifolds on Sgip, which defines the Melnikov function, is given by the zero order term of 
l/e[Ho{X''{9o,e)) - Ho{X'{9o, e))] |2|. In fact, if ^o, denotes the time evolution 

under H such that X^^'^{9o; 9q, e) = X'^/"(6'o, e), and Xo(t) denotes the time evolution under 
Hq such that Xo(0) = Xq, 

-±[X^I-{t-9Q,e)\ dt 

= -^[Xo{t-9o),t]dt + 0ie'). (12) 
Thus in the homoclinic case the Melnikov function can be written as 

Mi9o) = 11^^ [Xoit),t + 9o]dt, (^-^^ = {Ho,Hi}^ , (13) 

where { . , . } are the usual Poisson brackets. 

The implicit function theorem allows us to conclude that if M{9q) has simple zeros, 
then, for sufficiently small e the invariant manifolds intersect transversely for some 9q. On 
the other hand, if M{9q) is bounded away from zero, then the invariant manifolds do not 
intersect for all 6*0. 

Now let us take the map defined by system (0)-(@) on a arbitrary section. After a scale 
change in the p coordinate it reads 

Uk+i = Uk- Cul [pk + 0(4)] , (14) 
Pk+i=Pk-Cul[uk + {?>)], (15) 

where C = \^7i/Q. Here the standard Melnikov method breaks down because of the 
degeneracy of the saddle point. McGehee |24|, Xia [^, Dankowicz and Holmes [^] studied 



systems of this class in the context of the three body problem, where they established the 



fundamental results needed to support the Melnikov method: structural stability of the 
unstable periodic orbits (trivial in our case since the degenerated unstable periodic orbit 
remains fixed); existence of local stable and unstable analytic manifolds close to those 
of the unperturbed case solutions on the perturbed and the unperturbed manifolds 

approach to the unstable periodic orbit at a similar rate ||2^, Following the proofs step- 
by-step we can observe that all the essential hypothesis involved to achieve their results are 
also satisfied by the above system. Thus these statements apply to Eqs. ([T^ ) - ([T5|) allowing 



the expansion in Eq. (12) that justifies the use of Melnikov method in the present problem. 

Another important result whose standard form assumes hyperbolicity is the Smale- 
Birkhoff homoclinic theorem. It was given a formulation of this theorem that is valid for 



the degenerated problem of Sitnikov |26|, which is grounded on a suitable approximation 



of the linearization of the map in the neighborhood of the saddle point. Since such ap- 



proximation results in the same expressions for the case of Eqs. (p!^)-([T5|), we conclude 
that the Smale-Birkhoff homoclinic theorem applies to the above system. Thus transversal 
homoclinic intersections in our problem lead to the Smale horseshoe. 

Now we shall apply Melnikov method to Eq. (|l|) subjected to the perturbations (10). In 
the homoclinic case (0 < /3 < 1/16) we find for the Melnikov function 

M{9o) = / nr-''-^ cos[n{t + eo)]—dt (16) 

= -2nsm{neo)K{n), 

K{n) = / v''-^sm[nt{v)]dv, (17) 
Jo 

where we have fixed Xq = [v{0),p{0)] as the symmetry point on the loop and the integrand 
t{v) means the positive branch of Eqs. (|S|)-(§|)- Thus the Melnikov function has simple zeros 
as long as K{Q) ^ 0. With the change t ^ v we pass from an infinite interval in Eq. (|T6[) 
to a finite one in Eq. (|T^, and it allows us to study K{Q) using graphics. Although u 
(r = I/m^) is important to justify the Melnikov method, we have great freedom in the choice 
of a coordinate to study the results. The most simple one that is adequate to this end is the 
coordinate v = 1/r. 

The integrand of K{Q) is formed by the product of an oscillating function and a polyno- 
mial. Near the origin this oscillation is a rapid one, since t — > oo as the particle goes to the 
unstable periodic orbit. In Fig. 4 we show a graph of sin[f2t(f)] for (3 = 1/25, and several 
values of Q. For > 4 we will have more zeros in the interval shown in the figure. For f2 < 1 
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the curve will look like the one for Q = 1. Since the area under the curves are clearly not 
null the integral of sin[f2t(f)] is nonzero for j3 = 1/25 and < Q < A. The cases of interest 
are n > 2, where we will have a more favorable situation. We have transverse homoclinic 
orbits in all these cases. To better understand this behavior we show in Fig. 5 a graph of 
the integrand of K{Q) for f] = 1/25, = 3 and different values of n. For < f3 < 1/16 the 
graphics of sin[f2t(f)] will look like the one for f3 = 1/25, with almost the same upper bound 
for Q. See for instance Fig. 6, where we plot the graphics of sin[fit(f)] for the monopolar 
attraction center (/5 = 0) for the same values of Q employed in Fig. 4. Therefore we have 
transverse homoclinic orbits for all f3 limited by < /3 < 1/16. 

For these range of parameters the Smale-Birkhoff homoclinic theorem implies the exis- 
tence of a hyperbolic invariant set for which the action of an A^th iterate of the map has 
a symbolic dynamic equivalent to that of the Smale horseshoe map. Some important con- 
sequences of this result are sensitive dependence on initial condition (a characteristic of 
chaos); nonexistence of real analytic integral of motion (nonintegrable system); existence of 
infinitely many periodic orbits with arbitrary large periods (whose number increases expo- 
nentially with the period); capture of orbits by the system (in both directions of time). 

It is illustrative to see how the chaotic orbits look in the original noncompactified coor- 
dinates r,p. Due to the lack of an extra integral of motion the particle can have a highly 
erratic motion and have access to a two dimensional region of the phase space. The sensitive 
dependence on initial condition implies that the evolution of two infinitesimally near points 
in the space r,p can result in two completely different bounded orbits, in two completely 
different unbounded orbits, or even in one bounded and one unbounded orbit. Orbits that 
are bounded for all t < can go to infinite for t — > +oo, and orbits of particles coming 
from infinite can remain bounded for all t > 0. However, regular orbits are also present and 
in particular there is a family of periodic orbits. 

Let us consider the heteroclinic case (3 = 1/16. The Melnikov method obtained from 
Eq. ([T^) applies to each branch of the loop. The distance between X^^'^(±oo;9o,£) and 
X^/"(±oo;0o,O) is of order 0{e). Therefore Ho{X'/''{±oo;9o,e)) is of order ©(e^). Then, 
this term can be neglected and the Melnikov function reads 

M{9o) = 2nsm{neo)Ksin) - 2ncos(fi^o)^c(fi), (18) 

Ksin)= / v''-^sm[nt{v)]dv, K^{n)= / cos[^]t(^;)]rf^;. (19) 

Jo Jo 
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The orbits in the heterochnic loop are less symmetrical than in the homoclinic one and, 
consequently, the Melnikov function has simple zeros when at least one of the integrals of 
the previous formula is different from zero. The integrands of Kg and K^. oscillate very 
rapidly near the unstable periodic orbits. To better understand this behavior we show in 
Fig. 7 a graphic of the positive branch of Eq. (|^). The function t{v) has small values for 
a large range of values of f , and hence will result nonzero values for small Vt. Indeed, 
we show in Fig. 8 a graph of cos[nt(f )] for different values of VL. Since cosine is an even 
function, the same figure is valid for the negative branch. For > 3 we will have more zeros 
in the interval 1 < f < 4. For ^7 < 0.5 the curves will look like the one for = 0.5. For 
small values of n the integral Kc will be nonzero for < < 3. Due to the change of sign 
of cos[nt(f )] near v = 4, the upper bound for Q decrease with n. But it is clear from Fig. 8 
that for each n will exist a upper bound Qo{n) > such that Kc will be nonzero, leading to 
the presence of transversal heteroclinic orbits, for < < Qo{n)- 

IV. POINCARE SECTIONS METHOD 

The system (P) has four parameters, (3, n, Q, and e. In opposition to the Melnikov 
method, Poincare sections method is able to predict results only for fixed values of these 
parameters. However, the Poincare method can locate the regular and irregular regions and 
provide a qualitative idea of chaotic behavior. 

The perturbed phase space is in x and the maps defined on different values of the 
angular variable 6 = t mod 2-71 /Q are topologically conjugated. So we have restricted the 
study to sections built on = 0, but for a large number of different values of the parameters 
(3, n, Q, and e. We included Q values which would require a more elaborate numerical 
computation of the Melnikov function. In Fig. 9 we show the Poincare sections for some 
select values of /3, n, fl, and e. The general aspect of the sections is represented in this 
figure, with the predominance of irregular behavior near the destroyed invariant loop and 
the presence of a regular region near the center. 

A careful analysis of the Poincare sections reveals that chaos increases with n, which is 
natural since the major contribution of the perturbation comes from f > 1. In Figs. 9(a) 
and 9(b) we show the Poincare sections for n = 2 and n = 5, respectively, where (3 = 1/25, 
= 1, and e = 0.003. As a function of Q, the most chaotic behavior occurs for frequencies 
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of the same order of the angular frequencies of the orbits in the local minimum of Vejf, from 
(4/3)^ for (3 = 1/16 to 1 for j3 = 0. It is a reasonable result since for small fl the system 
is almost autonomous, and for large frequencies the particle feels only an average of the 
multipolar motion that goes to zero for Q — > oo. In Figs. 9(a), 9(c), and 9(d) we show 
Poincare sections for = 1, i7 = 0.1, and Q = 5, respectively, where (3 = 1/25, n = 2, and 
e = 0.003. Also, the chaos increases strongly with the "size" of the perturbation e. In Figs. 
9(a) and 9(e) we show the Poincare sections for e = 0.003 and e = 0.01, respectively, where 
P = 1/25, n = 2, and ^] = 1. 

The form of the orbits change with (3, since there is a change in the destroyed loop, but 
the relative chaotic area of the section is almost independent of the value of (3. In Figs. 9(a), 
9(f), and 9(g) we show the Poincare sections for (3 = 1/25, P = and (3 = 1/16, respectively, 
where n = 2, Q = 1 and e = 0.003. In particular, the case (3 = 1/16 presents two unstable 
periodic orbits and no higher chaotic behavior seems to be associated to it. 

V. FINAL REMARKS 

Making use of a compactification procedure and extensions of the standard Melnikov 
method and Smale-Birkhoff homoclinic theorem, we have detected transverse homoclinic and 
heteroclinic orbits and Smale horseshoe in periodically perturbed monopole + quadrupole- 
like potential. Moreover the compactified coordinates have proved useful to work out the 
Poincare section method. 

Although we have studied only one particular system, we stress that the approach is 
general and can be applied to a large class of open systems. The mathematical results depend 
essentially on the presence of the monopolar term in the unperturbed equations. Thus 
the same approach and with the same compactification applies to every inner-multipolar 
expansion model with nonzero monopolar contribution. In such cases we can allow non- 
Hamiltonian perturbations of the form f{u,t,e) in the equation corresponding to Eq. (^, 
where / is an analytic function of its variables, periodic in t, satisfying f{u, t, 0) = and of 
order 0{u^). 

These systems represent very general situations. In the gravitational case, for instance, 
they can model a potential due to a mass distribution moving periodically with reflection 
symmetry. This includes all planar potential written as Fourier expansion in time together 
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with an inner-multipolar expansion in space variable. A specific example is provided by 
a Newtonian binary system perturbed by gravitational radiation. This system models the 
long-term dynamical evolution of binary systems of stars due to the emission and absorption 
of gravitational radiation 

In electromagnetic systems we can have a different situation. The possible absence of 
monopolar contribution creates a difficulty in transforming the map equations into the ana- 
lytic form ([l^)-(|T5D. The solution may be found using the Cassasayas-Fontich-Nunes's result, 
which establishes the Melnikov method for systems with parabolic degenerated saddle points 



31| . A suitable compactification to treat the electromagnetic case is under study. 
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FIG. 1: The effective potential [Veff{r)] for /? = 1/10 (bottom curve), /? = 1/12, /? = 1/14, 
/?= 1/16, /3 = 1/18, and /? = (top curve), where the last one represents a monopolar potential. 
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FIG. 2: The level curves of the Hamiltonian Hq for /3 = 1/18. The homoclinic loop associated to 
{u,p) = (0,0) is the curve that contains that point. 
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FIG. 3: The level curves of the Hamiltonian Hq for /3 = 1/16. The curves that represent the 
invariant manifolds, which define the heteroclinic loop, associated to the points {u,p) = (0,0) and 
{u,p) = (2, 0) are the curves that contains those points. 



IF 



M-l 



-1 1 




0.6 1.5 2.4 



V 

FIG. 4: / = sin[J^t(t;)] for P = 1/25, and J7 = 1 (bottom curve), = 2, = 3, and = 4 (top 
curve). 
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FIG. 5: The integrand k = v''^'^ sm[nt{v)] of K{n) for /3 = 1/25, Q = 3, and n = 2 (bottom 
curve), n = 3 and n = 4 (top curve). 
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FIG. 6: / = sin[J7t(v)] for the monopolar potential (/3 = 0) with = 1 (bottom curve), = 2, 
= 3, and = 4 (top curve). 
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FIG. 7: The positive branch of Eq. (8): t{v) goes to ±00 dX v = and v = A, but has small values 
at almost every point of the interval < v < A. 
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FIG. 8: g = cos[^^^(^;)] for /3 = 1/16, and n = 0.5 (top curve), = 1, = 2, and = 3 (bottom 
curve). 
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FIG. 9: The Poincare sections for (a) /? = 1/25, n = 2, = 1, and £ = 0.003; (b) /? = 1/25, n = 5, 
O = 1, and £ = 0.003; (c) /? = 1/25, n = 2, O = 0.1, and £ = 0.003; (d) /? = 1/25, n = 2, O = 5, 
and £ = 0.003; (e) /3 = 1/25, n = 2, J7 = 1, and £ = 0.01; (f) /3 = 0, n = 2, Q = 1, and £ = 0.003; 
(g) /? = 1/16, n = 2, = 1, and £ = 0.003. 
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FIG. 9: (Continued) 
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